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Abstract. We show how to use the method of orthogonal poly- 
nomials for integrating, in the planar approximation, the partition 
function of one-matrix models with a potential with even or odd 
vertices, or any combination of them. 
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1. Introduction 



The method of orthogonal polynomials is a powerful technique for 
the non perturbative integration of matrix models over one or more 
matrices [0 in particular with even potential, i.e. with vertices with an 
even number of legs. Indeed, with even potential, the calculation sim- 
Q\ ', plifies both because the integrals are well defined and, as we shall see, 

the number of equations needed to solve the problem is smaller. On 
the other hand the model with odd vertices, in particular with cubic 
Ph. vertices is more natural in a number of problems; e.g. in the dynamical 

triangulation model of quantum gravity, where the random surface is 
given by a polyhedron with triangular faces, the order of the vertices 
appearing in the dual graphs is always three. Brezin et al. || solved the 
problem with cubic vertices using the saddle point technique. Bessis [Q 
introduced an alternative method (the orthogonal polynomial method) 
which to some extent appears more powerful e.g. in dealing with ma- 
trix model with more than one matrix variable [0. In particular, the 
orthogonal polynomial method has been proved useful in the treatment 
of a cubic vertex two-matrix model in the context of the Ising model 
on a random planar lattice. 

The purpose of this paper is to show, in a systematic way, how to 
extend the orthogonal polynomial method to arbitrary vertices, both 
even and odd and any combination of them. We shall follow the article 
of Bessis et al. generalizing some aspects to the case of odd vertices, 
in particular we shall recover, for the simplest case of cubic vertices, 
the result of M for spherical topology. Hopefully such a treatment can 
be extended to higher genus. 

The use of mixed vertices e.g. cubic plus quartic vertex, allows us to 
write a well defined i.e. convergent, partition function by adding to the 
cubic interaction a quartic term which makes the action bounded from 
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below, and thus the integral giving the partition function well defined. 
At the end one can take the limit when the coupling constant of the 
quartic vertex goes to zero. 

We start from the partition function 

Z N (g) = J dMe- trS < M > (1) 

where the integration is over an hermitian matrix of order N and where 
the action is given in general by 

- S(M) = ~M 2 + V(M) (2) 

k 

V(M) = V 9 [_ M j with k even and Re{g k ) < . (3) 

In the following we shall use also the variables = — f 2 — . As we can 

see the potential is composed of a combination of vertices, however, 
in order it to be bounded from above, and thus the integral be well 
defined, it is needed that the vertex with highest valence be even, and 
that the real part of the related coupling constant be negative. We are 
interested in the functions e^(g) 

ln^44 = TN 2 - 2h e h (g) (4) 

where 2 — 2h is the Euler characteristic of the oriented ribbon graphs 
to be summed in the perturbative expansion of the functions e/ l (g). 
Indeed, denoting such graphs with capital letters, each function 
admits the following expansion |I| 

o(A(G)) 



G connected of 
genus h 

where u»(Gr) is the number of vertices with valence % of the ribbon 
graph G and o(A(G)) is the order of the group of automorphisms of 
the graph []. Figures 0, ^| show the simplest connected ribbon graphs in 
the cubic and quartic case. Models with odd vertices will be regularized 
as explained above i. e. by adding a regulating even interaction which 
at the end is put to zero. We shall see explicitly how this method works 
for the cubic vertex. 



2 An automorphism of an oriented ribbon graph is defined in the following way. 
First let us identify the oriented ribbon graph as a common graph plus a cyclic 
ordering on the sets of half-edges attached to each vertex and then define an au- 
tomorphism of the oriented ribbon graph as an automorphism of the graph which 
leaves the ordering of each vertex unchanged. It's clear that the automorphism 
must send each vertex into a vertex of equal valence. 
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Figure 1. Second and fourth order connected graphs 
with cubic vertex. 
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Figure 2. First and second order connected graphs 
with quartic vertex. 
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2. The method of orthogonal polynomials 

A change of integration variables in (|I|) leads us to the integration 
over the eigenvalues A, of the diagonal matrix A 

Z N (g) = J dMe- trS ( M > = k H J fJdA, A 2 (A)e-£<S(*i) (6) 

where A(A) = HaK^'V ~ * s ^ ne Vandermonde determinant. We 
obtain the value of the constant ku using the results in ||: ku — 



1? jrpj. As we see the argument of the integral is the product of 

the Vandermonde determinant squared and a factorizable function of 
the eigenvalues, this feature makes the orthogonal polynomial method 
applicable. Let us introduce the measure dfx(X) = dAe" 5(A) , and the 
orthogonal polynomials P n (\) 

/+oo 
d f i,(X)P n (X)P m (X) = h n 6 rm (7) 
•oo 

where P n (X) is normalized by the condition that the coefficient of the 
term with highest degree equals 1 

P n {X) = \ n + . . . . (8) 

The polynomials P n (X) can be obtained in a constructing way e.g. by 
the Gram-Schmidt orthogonalization procedure from the monomials 1, 
A, A 2 , ... . A simple analysis of this procedure shows that the polyno- 
mials Pj have the well defined parity (— 1) J if the action S(X) is even. 
Every polynomial of degree n can be rewritten as a linear combina- 
tion of P m with m < n. The Vandermonde determinant in (^j) can be 
rewritten as 

JV 

A = det || Af 1 ||= det || P^X,) ||= J] i^j-iM (9) 

a i 

where the second equality is due to the fact that adding to a column a 
linear combination of the other columns does not change the determi- 
nant of the matrix; (— stands for the sign of the permutation a. 
We can take advantage of the coupling of the orthogonal polynomials 
due A 2 in @ to obtain the partition function in terms of the norm of 
the orthogonal polynomials 

N 

z N (s) = k H J2 (-i) p(CTl) (-i) p(CT2) n j MXi) P^-xiXdP^-i^i) 
= k H J2 (-l) p{ai) (-1) P(CT2) <W 2 II 'MO-i = k » m h ^ ■ ■ ■ h »~i 

a\, (72 i 

(10) 
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Let us rewrite this expression in a different form. The following equa- 
tion is valid 

AP n (A) = P n+1 (A) + A n P(X) + R n P n ^(X) (11) 

where the terms with index less than n — 1 are absent because after 
moltiplication by A they do not reach the degree n and thus are orthog- 
onal to XP n . For parity reasons when 5(A) is even A n vanishes. We 
shall refer to the preceding equation as the step equation because its 
repeated application enables us to calculate A*P n (A) using an analogy 
with all possible staircases i steps long. This method will be developed 
in the following section. Since 

h n +i = dfi(X) P n+1 XP n (X) 

J (12) 

= Jdfi(X) [P n+2 (A) +A n+1 P n+1 (A) +P n+1 P n (A)]P n (A) = R n+ ih n 

the partition function can be rewritten as 

Z N (g) = k H N\ h^R?- 1 . . . R 2 N _ 2 R N -! (13) 

where ho = J dAe -5 ^. Before passing to the limit for large N, we 
must compute 

hois) (14) 



1 Zjv(g) 1 v^/. n\ R n (g) 1 
En (g) = W ln ZMO) = N E l 1 - N ) ln lUO) + N ln 

v / n=l v 7 



h (0) 



We can prove that the last term on the r.h.s. of (|M| ) is negligible using 
the same perturbative method used for the partition function. The 
main difference is that h (g) is integrated over a scalar variable whereas 
partition functions such as Zjv(g) are integrated over a matrix variable. 
Expanding the action in ho and recalling the proof of the topological 
expansion JIJ, we notice the absence of the typical contribution due to 
the propagator deltas, namely a factor N for each face. Finally we have 
to add a factor N~ F which yields 

ln ^(o) " 2^ o(A(G)) • w 

uv ' G connected v v " 

Since X {G) -F = - E;= 3 - 1) < -1, we have 

which vanishes for large N. Thus 

e„(g) = >- E N(S ) = Jim i ± (l - £) in |M (17) 

n=l 
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3. The number of staircases 
We shall need, in the following, the quantities (3 l n defined by 

Kfa = J d/x(A)P n (A)A l P„_ 1 (A) . (18) 

We devote the present section to the calculation of the above integral. 
To compute \ % P n -\ we take advantage of an analogy with all staircases 
of i steps; where each step can go up, come down, or stay at the 
same level. The analogy comes from a repeated application of the step 
equation. After the integration only the staircases which end one step 
up, contribute. Each of them represents a product of factors: if a step 
is down from level n to the level n-lwe add a factor R n , and if it stays 
at the same level n we add a factor A n . Figure |3| shows an example of 
this kind of calculation. Since every coefficient A,-, Rj, is a function 



Pi = i 




= R „ +1 + R n + R n -,+ A 2 n + A 2 n ,+ A^ 

Figure 3. (3 l n computed from the staircases. 

of the index j it would be difficult to hand the final expression for j3 l n ] 
luckily, as we shall see, the planar limit ( N — > oo ) will enable us 
to neglect the differences among these quantities relative to different 
levels. In this limit we must compute the expression for (5 % n supposing 
that each step down yields a factor R, and each step that stays at the 
same level yields a factor A. Thus the question is: How many are the 
staircases of i steps whose final effect is to go up one step? Let j be 
the steps of type A, then the other i — j are divided in p steps up and 
p — 1 steps down so that i — j + 2p — 1. Without the A steps there 
are (^p -1 ) staircases of 2p — 1 steps whose final effect is to go up one 
step. Inside these staircases we want to insert the remaining j levels of 
type A: there are 2p places where they can be inserted, and, for a fixed 
staircase, there are ^ 2p+ ^~ ) choices. Finally the number of staircases 
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of i steps whose final effect is to go up one step is 



|J (i-2p+l)!p!(p-l)! ' (19) 

where [] stands for the integer part, and, denoting by (3 l the continuum 
value of (3 l n) we have 



where the tilde reminds the replacement Aj — > A, Rj — > R. The values 
of (3 l for the first few % are 

p 2 = 2A 

/3 3 = 3A 2 + 3R (21) 
/3 4 = AA 3 + Y1AR . 

Analogously we define 

h n+1 y n = J dfi(X) P n+1 (X)\ l P n ^(X) . (22) 



By the same technique used for (5 l n we find 

" =2 \i-2 P + 2)\ P \( P -2)\ 
Finally we define 



i+li 



* - ^-'=E (i _ 2p ( ;; ) ;^_ 1) ,^ - w g- 1 • 

(24) 

The values of 5* for the first few i are 

5 3 = A 2 + 2i? 

5 4 = A 3 + 6AR (25) 

5 5 = A 4 + Y2A 2 R + 6R 2 . 



4. Derivation of the continuum equations 



In this section we shall examine the continuum limit N — > 00, which 
will allow us to write a simple expression for the generating function 
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e 2 (g) of the planar graphs. This will also justify the replacement Aj — > 
A, Rj — > R used in the previous section. Let us consider the identity 

nh n = J d/i(A) \P' n {\)P n {\) 

= J d/i(A) P' n {\) [P n+ i(A) + RnPn-l(X) + A n P n {\)} 

= R n j d\e- s ^P' n {\)P n ^{\) (26) 

= R n J dXe^S'(X)P n (X)P n ^(X) 

k 

= (l-Y,9i?n 1 )h n Rn, 



i=3 



where in the last but one equality we have integrated by parts and 
in the last equality we have used the definition of j3 l n . Thus we have 
obtained the first recursion relation 



n 

i=3 



From this equation we infer in particular that: R n (0) = n. We want to 
find a second recursion relation which relates the coefficients A n and 
R n . We observe that: 

J d\e- s ^\P n (\)P: +1 (X) = J dXe~ s ^P n (X)XS'(X)P n+1 (X) 

k (28) 

= (A n + A n+l - ^ 9iPn+l) K+i ■ 
i=3 

But 

J d\e- s ^\P n (\)P; i+1 (\) = nA n h n + R n Jd\e- s ^P n ^{\)P ! n+1 {\) 

= nA n h n + R n J dAe" s W 'P n _ 1 (X)S'(X)P n+1 (X) 
k 

= nA n h n - h n+1 Rn ^ gcf' 1 . (29) 

i=3 

As a result, the second recursion relation is 

k k 
(A n + A n+1 - 9iP l n+i) Rn+i = nA n - R n+1 R n ^ gcf' 1 . (30) 

i=3 i=3 
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Now, we extract the planar case taking the limit iV — > oo. Let us 
introduce the substitutions 



n 
N 

~w 

A n 



x (31) 
R(x) (32) 
A(x) (33) 

f N 

to obtain, taking into account the power of N contained in g^ the two 
continuum equations 

k 

x = R(x)(l-^ gi ?- x (xj) (34) 

k 

A(x) = Y^9rS\x), (35) 

i=3 

where j3 l {x) and 5 l {x) are expressed in terms of A(x) and R{x) as given 
by eqs. ( pOj , [24]) . One finds, from eqs. (|27|, ^0|), that the continuous 
solution A(x), R(x), is related to the coefficients A n and R n by 

A 

AT 

We are interested only in the continuous solution, indeed, recalling that 
_R n (0) = n, we have 

#n(g) i2n(g)/nN"l 



__*(£) +0(jn 



a — - ( 36 > 



1 



(37) 



and the function eo(g) can be rewritten, in the limit N — > oo, as 



e (g)= / * da: (1 - s) In I . (38) 



o 



X 



We have therefore reduced our problem to the one of obtaining R(x) 
(or A(x)) from the system of continuum equations (|34], |35]), and to 
perform the integral in eq. (|38|). 



5. The solution for the cubic vertex 

In order to understand the solution for the cubic vertex it is useful 
to recall the main features of the even vertex case. For even vertex 
A(x) = and thus the first continuum equation for R(x), 

- = R (l - E ( p ) RP ^) ( 39 ) 

p— 2 
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suffices. The quartic case can be explicitly integrated to obtain 

1 1 00 (2k -j \ | 

eoM = - In o + — (a - l)(o - 9) = $>0 4 ) fc fc \ ~ ^ (40) 

fe=i ' ^ 

with 



WT^ = y. t (2t-l)l 

Recalling the formula for the topological expansion fl5|), one has the 
interesting equation 

V 1 = t ^k-l)\ (42) 

^ o(A(G)) A;! (A; + 2)! 1 ; 

G planar, connected, 
with k quartic vertices 

that can be checked for k — 1 and k = 2 using the contents of figure 0. 
The radius of convergence is 1/12 and g4 C = 1/12 is the critical point. 
For #4 — > g 4c one obtains the critical behavior 
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eo^) ~ (94c - 94) 2 , (43) 

which is particularly interesting for 2D-Gravity. 

Coming now to the cubic vertex case, in order to cure the lower 
unboundedness of the action, one can add a quartic term with negative 
#4, with the idea to take eventually the limit g± — > 0. On the other 
hand it is clear that we can just take g 4 = in all the expressions which 
turn out to be well defined in that limit. We shall see that the integral 
( pIHD is well defined for R(x) given by the solution of the continuum 
equations Q3"4], 1351) with g 4 = 



- = 1 - 2g 3 A 

R (44) 
2g 3 R = A-g 3 A 2 . 

In fact, let us introduce the new variable a = —g 3 A related to x by 

2g 2 3 x + a(l + tr)(l + 2a) = . (45) 



The function a(x) = a(g 3 x) is the solution of ( [45] ) which vanishes for 
x = 0; indeed when g 3 = the potential has no longer odd vertices 
and then A(x) = =>- er(0) = 0. Our function eo(gs) may be rewritten, 
going over from the variable x to the variable a and integrating by 
parts 



(46) 



2 v " 3(l + <7i)(l + 2<7i; 
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where o\ = cr(gf ) is the solution of 

2 ^ + (T (l + a)(l + 2a) = 0, (47) 

which vanishes in g 3 = 0. <i\ can be expressed as an expansion in 
powers of g% using Lagrange theorem, obtaining 

ai --ih i!h:} r(t + i)r(l(t + i))- (48) 

Except for some factors, due to different definitions, our results coincide 
with those of Brezin et al. ||. The power expansion series for the planar 
generating function is 

and, recalling the topological expansion for eo(<?3), we reach the formula 

V 1 = 8 *i r ( 3fc / 2 ) ( , 0) 

^ o(A(G)) 2T(H3) m/2 + 1)' y ' 

G planar, connected, 
with 2k cubic vertices 

Such formula can be checked for k — 1 and k = 2 using the contents of 
figure 0. The asymptotic behavior of the expansion coefficient is 

,1 r(3fc/2) e 3 , r-.u, v , N 

8 - -T-, ^-77^ r =(12V3) k ks . 51 

2 r A; + 2) r (k 2 + 1) ; K J 



The radius of convergence of the series is 1/ \^12~\/3 and g? jC = 

l/\/l2\/3 is the critical point. From the asymptotic behavior 

we conclude that the critical exponent remains unchanged from the 

quartic case. 

In 2D-Gravity, where the continuum surfaces are replaced by poligo- 
nalizations, such a result is a check of the independence of the partition 
function, in the limit of infinite number of vertices, of the kind of polig- 
onalization one chooses to approximate the continuum surfaces [[Tfl. 

6. Conclusions 

In dealing with matrix models usually one encounters matrix models 
with even potential so the question naturally arises if there is some 
obstruction to the odd vertex case. In this paper we have shown that, 
even if, in the odd vertex case, the original partition function is ill 
defined, the method of orthogonal polynomials can be often applied in 
its most naive form, that is ignoring all convergence problems. This is 
justified by adding a regulating even vertex to the odd one, and taking 
eventually the limit for its coupling constant going to zero. We have 
extended the orthogonal polynomial method to any combination of odd 
and even vertices, writing the two needed continuum equation. The 
explicit application to the cubic vertex case has been given, recovering 
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the result of Brezin et al. |J . An explicit integration of 3+4 or 5 vertex 
case appears feasible along these lines and would be a useful check of 
the universality of the critical behavior. 

The general setting explained here can be readily developed, in the 
planar case, also for two-matrix models with coupling in the form of 
the Itzykson-Zuber formula || 0, the cubic case being already solved 
in ||. Further extension can be developed in higher genus cases e.g. 
in the cubic case for the torus topology. 
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